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Abstract 

Electron paramagnetic resonance (EPR) measurements on single crystals of 
the molecular magnets Feg and Mni2 reveal complex nonlinearities in the 
linewidths as functions of energy eigenstate, frequency, and temperature. Us- 
ing a density-matrix equation with distributions of the uniaxial anisotropy pa- 
rameter D, the Lande g factor, and dipolar interactions, we obtain linewidths 
in good agreement with experiments. Our study shows that the distribution 
in D is common to the examined molecular magnets Feg and Mni2 regardless 
of the qualities of the samples. This could provide the basis for a proposed 
tunneling mechanism due to lattice defects. The distribution in g is also quite 
significant for Mni2. 
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I. INTRODUCTION 



Molecular magnets such as Mni2^ and Fcg^ have recently drawn vigorous attention 
because of the macroscopic quantum tunneling (MQT) of their magnetizations at low 
temperatures^'^ and their possible applications to quantum computing.^ These materials 
consist of many identical clusters with the same magnetic properties and characteristic en- 
ergies. Each cluster is made up of many different species of ions and atoms, with a total 
spin angular momentum in the ground state of 5'=10. The clusters have strong crystal-field 
anisotropy and, thus, a well-defined easy axis, and the magnetic interaction between different 
clusters is weak. Many competing models have been proposed to explain MQT in molecular 
magnets: higher-order transverse anisotropy,^ thermally- assisted quantum tunneling,^'® the 
Landau-Zener effect, and dipolar interactions with dynamic hyperfine fields. So far, the 
MQT in these materials has not been completely understood, and new experimental and 
theoretical approaches are needed. 

In this paper, we present both previously unpublished experimental data and a theo- 
retical model for the electron paramagnetic resonance (EPR) linewidths of single crystals 
of the molecular magnets Fes and Mni2 when the field is applied along the easy axis. To 
understand the measured linewidths, we use density-matrix equations, assuming that the 
uniaxial anisotropy parameter D and the Lande g factor are randomly distributed around 
their mean values ("D-strain" and "^f-strain" effects^^) due to crystalline defects or impuri- 
ties in samples. By adjusting the standard deviations of D and and the strength of the 
dipolar interactions, these equations are used to obtain theoretical linewidths in excellent 
quantitative agreement with our experimental data. Our theoretical study also shows that 
for the examined Feg sample the distribution in D and the dipolar interactions are crucial 
to understand the linewidths, while for the Mni2 sample the distributions in D and g are 
more important than the dipolar interactions. These results imply that the distribution in D 
seems to be generic in the molecular magnets Feg and Mni2, although the standard deviation 
of D itself varies between samples. This point assures that inherent defects or impurities in 
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samples can be checked through the distribution in D by observing the hnewidths of EPR 
spectra, and forms a basis of a recently proposed dislocation-induced tunneling mechanism. 

The experimental details are presented in Sec. II. Section III describes some typical 
EPR spectra and observed features, while Sec. IV presents our theoretical model. Section 
V discusses the theoretical results in comparison with the experimental data. Section VI 
describes our conclusions. 

II. EXPERIMENTAL DETAILS 

All the measurements were made on single crystals. Mni2 was synthesized according 
to the procedure described by Lis.^ Fcg was synthesized using Wieghart's method.^ Single 
crystals of both materials were grown via slow evaporation. Typical crystal dimensions used 
were IxlxO.Smm^. EPR measurements were made for a magnetic field along or close to the 
easy axis using the multifrequency-resonator-based spectrometer described earlier. ^^^^^ The 
(transverse) excitation frequency was varied between 45 and 190 GHz. The temperature 
range was from 1.8 K to 40 K, with a precision of 0.05 K. The linewidths were measured as 
functions of Zeeman field, EPR frequency, and temperature, for crystals of different quality, 
over a wide enough range to probe energy-level broadening due to lattice defects and inter- 
cluster dipolar interactions. 

III. EXPERIMENTAL RESULTS 

Figure 1 shows a typical EPR spectrum of Fes for a field approximately along the easy 
axis at a resonance frequency of 133 GHz and at T=10 K.^^ Figure 2 shows a typical EPR 
spectrum of Mni2 for a field along the easy axis at a resonance frequency of 169 GHz and 
T—25 K. The following interesting features are observed: the EPR linewidths are on the 
order of several hundred to a thousand gauss for Fes (one to a few thousand gauss for 
Mni2) [Figs. 1, 2, 3(a),(b),(c), and (d)]; the linewidths increase non-linearly as a function 
of the absolute value of the energy eigenstate Mg [Figs. 1, 2, 3(a),(b),(c), and (d)]; the 
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minimum linewidths are observed for transitions involving the Mg—O eigenstate, with a 
weak asymmetry in the widths about Ms=0 [Figs. 3(a) and (b)]; and resonances involving 
the same energy eigenstates become narrower upon increasing the measurement frequency 
or field (compare Fig. 3(a) with 3(c)). 

IV. MODEL 

For clarity, the model for Feg is described separately from that for Mni2. 

A. Fes 

To obtain the theoretical hnewidths of the power absorption for Feg clusters, we consider 
a single-spin system with S=10 in a weak oscillating transverse magnetic field. We choose 
the easy axis to be the 2;-axis. Since the Feg clusters have an approximate D2 symmetry, 
the lowest-order ground-state single-spin Hamiltonian is^^ 

Ho = -DSl - E{Sl - Si) - g^iBH.S, , (1) 

where D=0.2S8kB {ks is Boltzmann constant) and the transverse crystal-field anisotropy 
parameter, is 0.043A;b.^^ Here Sa is the a-th component of the spin angular momentum 
operator, g is the Lande gi-factor which is close to 2, is the Bohr magneton, and is 
the longitudinal static applied magnetic field. In zero field, ignoring the small transverse 
anisotropy terms, the eigenstates [1/ \/2){\Ms) ± |— M^)) are degenerate for each M^, and the 
lowest energy states are (l/-\/2)(|10) ± | — 10)). These degeneracies are lifted by applying the 
field, H^. When ~ D{m -\- m')/gfiB, where m and m' are eigenvalues of in units of H 
(m, m'=— 10, • • • , 10), there exists a pair of quasi-degenerate energy levels, m and m'. When 
Hz is large enough, the eigenvalue of S^, Mg, is a good quantum number since the transverse 
terms are much smaller than the longitudinal terms. In such magnetic fields, the transverse 
terms, E{Sl — Sy), cause quantum tunneling between the two quasi-degenerate energy levels. 
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Next we introduce an interaction V{t) between the spin system and an oscillating transverse 
magnetic field with angular frequency a;=27r/: 

V(t) = ^{e'^' + e-^') , Vo oc H,S, . (2) 

We treat V{t) as a small perturbation to Ho- 

Besides the interaction V{t), the spin system interacts with the environment, such as the 
thermal fluctuations of the clusters. Through this spin-phonon coupling, the spin system 
can relax to re-establish thermal equilibrium. Since the times of interest are much longer 
than the correlation time of the phonon heat bath, the heat bath is considered always to 
be in thermal equilibrium. Thus we can integrate out all degrees of freedom of the bath to 
obtain an equation of motion for the density matrix p{t) of the spin system:^'^ 

(3) 

lm',m = , Wm ^ 2^ Wkm , (4) 

where the subscripts represent eigenstates of the longitudinal part of Ho, p{t)m'm 
= (m'|p(t)|m), and Wkm is the transition rate from the m-th to the /c-th eigenstate, which is 
determined by the spin-phonon coupling.^ 

We consider the case that the frequency u of the oscillating field is fixed while is varied 
to induce a resonance. With the selection rule AM5=±1, the power absorbed between the 
Ms and the Ms - I level is 

dS- 

where £ms is the energy of the Mg level, and [pM^.Mjrad is the change with time of the 
population in the level Ms caused by the first term in Eq. (3). We solve Eq. (3) for [pMs,Mjrad 
up to first order in Vq near resonance \uj ~ {^^Ms ~ ^m^-i)/^]- In this limit, among the 
many off-diagonal density-matrix elements, only p{t)Ms-i,Ms [= p{t)*Ms Ms-i\ contributes to 
[PMs.AfJrad- The matrix element p{t)Ms-i,Ms obtained by solving Eq. (3). Consequently, 



only two neighboring energy levels are involved for each resonance peak. An increase of 
the population in one of the two energy levels results in a decrease of the population in the 
other. So the power absorption becomes^° 



— V /^[J^z)[PMs,Ms-PMs-i,Ms-i) > 



dt 

1,2. 

l,M: 



V ^ \mv„\M. - 1)1. ^ ^-mM.-i) 

9Pb 

where A{Hz) is a Lorentzian lineshape function, ifj-es is the resonant field, and 
hlMs-i,Ms/ 9Pb gives a linewidth due to the finite lifetime of any excited state. As Mg 
decreases, the resonant field increases. The linewidth determined by Jms-i,Ms is on the 
order of several to several tens of gauss at temperatures below several tens of kelvin, and 
it increases as Mg decreases because of the corresponding shorter lifetimes. However, the 
experimentally observed linewidths decrease as Mg decreases until the Mg—1 and 
Mg—0 — > —1 transitions (hereafter denoted Mg—1 and Mg—0, respectively), and are on the 
order of several hundred to a thousand gauss in the same temperature range [Figs. 1 and 3 
(a),(b), and (c)]. 

To resolve these large discrepancies, we first assume that D and g are independent random 
variables with Gaussian distributions centered at 0.288A;b and 2.00, with standard deviations 
ao and ag, respectively. Such distributions can be caused by impurities or crystalline defects 
in the macroscopic sample, which result in different clusters seeing slightly different values 
of D and g. Then we calculate the average power absorption at a fixed frequency and 
r=10 K by averaging Eq. (6) over the Gaussian distributions using Mathematical^, to obtain 
the linewidth as a function of Mg. Our numerical calculations show that the distribution 
in D causes the hnewidths to decrease linearly with decreasing absolute value of 2Mg—l, 
with a slight rounding close to the linewidth minimum (Ms=l,0), where the D-strained 
linewidths approach the intrinsic lifetime broadening. On the other hand, the distribution 
in g makes the hnewidths increase with decreasing Mg because the resonant field increases 



6 



with decreasing Mg. These trends can be predicted from the expression for the resonant 
field Hj-es in Eq. (6) because the intrinsic hnewidths are substantially smaller than the 
measured linewidths (with the exception of the lowest absolute Ms transitions). However, 
the density-matrix equation (3) is needed to understand the temperature dependence of the 
linewidths. 

Next, we consider the effect of the inter-cluster dipolar interactions, in which each cluster 
experiences a net magnetic field comprised of the applied field and the dipolar fields from 
the surrounding clusters. At high temperatures {kBT^Tno) these dipolar fields fluctuate 
randomly about zero, giving rise to a dipolar contribution to the EPR linewidths. At low 
temperatures, the ground state becomes preferentially populated, leading to nonzero dipolar 
fields which are the same at the site of each cluster. This leads to a hne narrowing upon 
increasing i^res or the EPR frequency, or reducing the temperature. We use the multi-spin 
Hamiltonian that consists of TIq and V{t) for each cluster, and magnetic dipolar interactions 
between different clusters. We neglect the spin-phonon interaction because the spin-phonon 
relaxation time is very large compared to the spin-spin relaxation time. We truncate the 
original multi-spin Hamiltonian to leave only terms which commute with the z component of 
the total spin angular momentum Szj, where the sum runs over all clusters, because of the 
selection rule AM5=±1 where Mg is an eigenvalue of J2j Szj- Then the dipolar interaction 
is modified to n"^ = Eik>j Ajk{Sj ■ - 35,,-5,fe), where A^k = [{gfiB)V2r%]{3Cfk - 1)- Here 
rjk is the distance between clusters j and k, and Cjfe is the direction cosine of rjk relative to 
the 2;-axis. Assuming that J2j Vj{t)<H'^'^^Y.j ^oj, 

we neglect Y^j the energy of the 

multi-spin system and treat 7^^^ perturbatively. Following the formalism in the literature^^ 
for the 5'=10 system, we calculate to first order in Ti'^P the root-mean-square deviation of 
the resonant field from the center value i^res in Eq- (6), assuming that the line is symmetric 
about the center value. We do not account for the exact geometry of the sample for this 
calculation. We use a lattice sum, assuming that dipoles are distributed on a simple cubic 
lattice, and the magnetic anisotropy axis is not particularly aligned along any of the edges 
of the lattice. The specific angles of the magnetic anisotropy axis with the lattice axes do 

7 



not appreciably change the strength of the dipolar interactions. The dipolar interactions 
make the linewidths decrease as Mg decreases because the stronger resonant fields lead to a 
more polarized system. 

In this study, we emphasize that the line broadening effect due to hyperfine fields does 
not appear to be a significant factor. Indeed, for Fcg, there is no hyperfine splitting due to 
the metal ion, since the dominant Fe isotope, ^^Fe, has no nuclear spin. 

The overall features of the hnewidth, as a function of M^, are determined by the competi- 
tion among these three effects: i?-strain; g^-strain; and dipolar interactions. We have varied 
ao, CFg-i and the effective distance between the neighboring dipoles, within experimentally 
acceptable ranges in order to fit the experimental data. 

B. Mni2 

Since the molecular magnet Mni2 has tetragonal symmetry, the ground-state single-spin 
Hamiltonian is, to lowest order, 

Tio = -DSl - CSl - gi^BH.S, (7) 

with D—0.55kB, C—l.VI x IO^^/cb, and g—l.QA?^ We perform the same analysis as in the 
previous subsection for the available Mni2 data (/=169 GHz, T—2b K, the field along the 
easy axis). See Fig. 3(d). The differences are as follows: the resonant field is modified to 

^ fiuj-D{:2M, - 1) + C(U/^ -6^7; + [M„ - 1) _ 

-"res = ) 

dipoles are assumed to be distributed on a bcc lattice; and the anisotropy axis is along the 
longer side of the sample. 

For Mni2, the hyperfine effect must be weak, because there is no resolved hyperfine 
splitting in the EPR spectra. We attribute this absence of hyperfine splitting to exchange 
narrowing due to fast electron spin exchange over the metal-ion framework. Otherwise one 
would expect a hnewidth of over 1000 G, since each ^^Mn nucleus (/=§) is expected to exhibit 
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a hyperfine splitting of about 200 G.^^ Thus 12 almost equivalent ^^Mn nuclei, together with 
the acetate protons alone should lead to an overall width of over 1000 G. Notice that the 
measured smallest width {Ms=l in Fig. 3(d)) is about 1000 G, which increases to over 2000 
G toward the transition Ms=10— >9. The hyperfine effect would not be expected to change 
significantly with the electronic spin quantum number M^. We thus consider the hyperfine 
effect as a lumped parameter representing the residual linewidth, independent of Mg. 

V. RESULTS 

For the Feg sample examined, the calculated linewidths agree well with the experimental 
data at the measured frequencies (/ = 68, 89, 109, 113, 133, and 141 GHz) at T=10 K, using 
ao^O-OiD and df^l2 A. Within acceptable ranges of an, o-g, and d, no other combination 
of parameter values is able to produce the same quality of fit simultaneously to all of the 
frequencies studied, although other parameter values may produce good fits for a particular 
frequency. Figures 3(a), (b), and (c) show the experimental and theoretical linewidths for 
Fes at three frequencies. Our calculations show that the £)-strain and the dipolar interactions 
are equally important for the linewidths of the sample, while the gi-strain does not contribute 
significantly. It was previously speculated that Mg dependence of the EPR hnewidths for Fcg 
and Fe4 was due to the £)-strain only.^^ However, this cannot be explained by the £)-strain 
alone, because the distribution in D only results in a linear Ms dependence of the linewidths. 
The possibility of the presence of dipolar fields was mentioned in a zero-field millimeter- wave 
experiment for Feg.^^ The weak asymmetry of the hnewidths about Mg—l and Mg—O, shown 
in Figs. 3(a) and (b), can be understood by combining the D-strain effect with the dipolar 
interactions. The D-strain produces linewidths symmetric about Ms—1 and Ms—0, while 
the dipolar interactions yield narrower linewidths with smaller Ms transitions, which are 
observed at higher resonant fields. For a given Ms transition, the linewidths decrease as 
the frequency increases because of the dipolar interactions. Compare Fig. 3(a) with 3(c). 
The distribution in D does not give frequency- or resonant-field dependent linewidths, but 
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the dipolar interactions do. This is another reason that the distribution in D alone cannot 
explain the measured linewidths. We also considered distributions for other higher-order 
crystal-field anisotropy parameters {E or fourth-order coefficients) and found that they do 
not contribute substantially to the linewidths. 

For the Mni2 sample examined, the calculated linewidths [Fig. 3(d)] are in a good agree- 
ment with the experimental data for ag^O.OQSg, and (i~14 A. Figure 3(e) shows 
separately the line broadening due to the D-strain, g^-strain, and dipolar interactions with 
the chosen standard deviations and the effective distance between dipoles. For the Mni2 
sample, the distribution of D is wider than that for the Fes sample. The D-strain effect is 
still dominant, but a g^-strain effect must be included as well to understand the linewidths 
for smaller Mg transitions. As shown in Fig. 3(e), the linewidths for and Mg—l^O 
of over 1000 G cannot be explained by the D-strain and the dipolar interactions alone. It is 
essential to invoke the g^-strain effect. The dipolar interactions do not play as significant a 
role as for the Fes sample. 

VI. CONCLUSIONS 

We have presented our measured EPR spectra for the single-crystal molecular magnets 
Fes and Mni2 as functions of energy eigenstate Mg and EPR frequency when the field is 
along the easy axis. We also have calculated theoretically the EPR linewidths to quan- 
titatively compare with the experimental data, using density-matrix equations along with 
Gaussian distributions in D and g and the dipolar interactions. Our calculated linewidths 
agree well with the experimental data. Our analysis shows that for the examined Fes sample 
the distribution in D and the dipolar interactions play key roles in explaining the linewidths, 
while for the Mni2 sample the distributions in D and g are more important than the dipolar 
interactions. From these, we draw the conclusion that the £)-strain is universally important 
for the EPR hnewidths in molecular magnets, although other effects, such as g'-strain and 
dipolar interactions, can also contribute, depending on the material. These results imply 
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that the assumption made in the tunnehng modeP^ may be correct. A distribution in D 
may also facihtate multi-photon transitions through virtual states which are necessary to 
implement a quantum algorithm for the molecular magnets Mni2 and Feg.^ Additionally we 
speculate that the £)-strain effect may contribute to the widths of resonant steps in magne- 
tization hysteresis loops. A more detailed report, including the temperature dependence of 
the linewidths at fixed frequency for both molecular magnets, will be published elsewhere. 
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FIG. 1. Normalized power absorption (in arbitrary units) for Feg as a function of magnetic 
field at r=10 K with the field almost along the easy axis {6 = 9°y^. Every other resonance peak 
is marked by the energy level from which the spin system is excited. For example, the leftmost 
peak is for the transition M5=10^9. 
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FIG. 2. Normalized power absorption (in arbitrary units) for Mni2 as a function of magnetic 
field at T=25 K with the field along the easy axis. Every other resonance peak is marked by the 
energy level from which the spin system is excited. The leftmost full resonance peak around 1 T 
is for the transition Ms=6— ^5. 
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FIG. 3. Experimental (filled circles) and theoretical (open circles) Full Width at Half Maximum 
(FWHM) of the resonance lines vs Mg for (a),(b),(c) different frequencies at T=10 K for Fcg, and 
(d) at r=25 K for Mni2. (e) line broadening due to -D-strain, g'-strain, and dipolar interactions 
vs Ms for the Mni2 sample in (d). The standard deviations of D and 5, and the effective distance 
between dipolcs for (e) are given below. Here For Fes (Mni2), the fit is best when the 

standard deviation of D is O.OII? (0.02D and the standard deviation of g is O.OO85) and the 
effective distance between dipoles is 12 A (14 A). 
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